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I We perform a systematic operator product expansion of the most 

I general form of the nucleon scattering tensor W^^i, including electro- 

I magnetic and weak interaction processes. Finite quark masses are taken 

I into account and a number of higher-twist corrections are included. In 

I this way we derive relations between the lowest moments of all 14 struc- 

I ' ture functions and matrix elements of local operators. Besides repro- 

^ I ducing well-known results, new sum rules for parity-violating polarized 

structure functions and new mass correction terms are obtained. 
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I. INTRODUCTION 



Deep inelastic lepton-nucleon scattering (DIS) provides the cleanest method to un- 
derstand the internal structure of the nucleon Presently such measurements 
performed at HERA, CERN, SLAG, and other accelerators concentrate on electro- 
magnetic DIS. Theoretically weak interaction processes can provide valuable addi- 
tional information [Q, allowing e.g. a flavor decomposition of structure functions. 
Such experiments are extremely difficult and require careful and detailed analysis, 
but the possible scientific progress justifies a serious consideration of them 

All DIS experiments aim at measuring the structure functions of the nucleon. In 
the framework of operator product expansion (OPE), moments of structure functions 
can be related to nucleon matrix elements of well defined local operators which have to 
be calculated nonperturbatively, e.g. in lattice QCD 0]. In some cases those matrix 
elements occur in other strong interaction processes and sum rules can be derived. 
The polarized Bjorken sum rule ^ , relating nucleonic /3-decay and the first moment 
of the polarized spin structure function / dxgi{x), is an especially interesting case. 
Such sum rules offer a convincing way of testing perturbative and non perturbative 
QCD. 

With the accuracy of QCD tests reaching the per cent or even per mille level, 
the higher- twist contributions to DIS become relevant. While quite a number of 
theoretical investigations on their form exists scattered in the literature, no system- 
atic presentation of all contributions up to the second order in the expansion of the 
fermion propagator within the formalism of OPE has been given so far. As compared 
to the interactionless quark parton model used in the analysis of @], the OPE for- 
malism is of general validity and includes the QCD interactions. This can be seen in 
the fact that within the framework of the quark parton model obtains for example 
g2{x) = 0, while the OPE predicts Jq dxx'^gl^x) being different from zero. 
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This paper gives a complete determination of the lowest moments of the 14 nu- 
cleon structure functions which parametrize the most general form of the nucleon 
scattering tensor (only restricted by Lorentz and time-reversal invariance). The pa- 
per is organized as follows: In section II we define all relevant quantities and follow 
the standard procedure of OPE to isolate contributions of well defined twist. This 
section is supplemented by appendix A in which the decomposition of all relevant op- 
erators is performed explicitly. Using these results we analyze 7 and Z° exchange in 
section III and exchange in section IV. Our final results are presented in section 
V and the conclusions in section VI. A discussion of the contributions from cat-ear 
diagrams is delegated to appendix B. 

Our main results are new higher-twist and quark mass corrections to the parity 
violating polarized structure functions. 



II. DEFINITIONS AND SPIN DECOMPOSITION 

DIS processes can be distinguished according to the exchanged boson, namely 7, 
or W^. The cross sections for these processes are 



dE'dn EM (Q2) 



dE'dn EMiQ-' + iMzoYY \ ^ MlJ\ ^ M2 ' ' ^ ^ 
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dE'dn EM{Q^ + {Mw±)T V M^^)[^ M^^J ' 

Here E is the energy of the incoming lepton, E' the energy of the outgoing lepton, 
and M is the nucleon mass. is the sohd angle of the outgoing lepton. As usual q is 
the photon, respectively boson momentum with = —q'^, a is the electromagnetic 
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coupling constant. The weak interaction coupling constants are absorbed into the 
lepton (L^iy) and nucleon (Wfj,,^) scattering tensors. 
These tensors can be written in the following form: 

a' 

\pX) . (4) 

The index Pr = 7, Z^, W~^, W~ indicates the process under consideration. The 
lepton current is labeled by L and the hadronic one by H. is the four-momentum 
of the incoming lepton and its spin, k'^ and a'^ are four-momentum and spin of 
the outgoing lepton. The nucleon four-momentum is p^, {p^ = M^) and A = ±1/2 
denotes its polarization (cf. figure 1). Let 5*^ be its spin vector , then according to 


Wl,:^\q,p,S) = E(A|p|A) {Wl,^^\q,p))^ , (5) 

A 

with the spin density matrix 

and = {0,S) in the rest frame of the target (S*^ = — M^, S ■ p = 0). We use 
covariant normalization for the lepton states {ka\k'a') = 2kQ{2nY6{k — k')6„y, and 
{pS\p'S') = 2po{2TTY5{p — ff)5s,s' for the nucleon states in Wj^^^^q^p, S). The lepton 
currents read 

e 

I 

e,v 

i^'^^'KO =T.mi,iV'- + A'-75)e(0 , (7) 
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d'^e^^HpM 



where / G {e ,/i ,r , z/g, z/^, z/^}, e G {e , /i ,r } and v G {z^e, t'/i, z^r}- The nucleon 
current is written as 



jI^'^^HO = 9«(07A.(K!i + Kdl^Md^i) + further flavor combinations , 



(0 = E^/(07M(^/ + ^/75)g/(0 



(0 = (IdiOilp^'^du + ^du75)gn(0 + further flavor combinations 



The constants V and A follow from the standard model (see for example [TT|]). Qf 
is the electromagnetic quark charge and / runs over all flavors. All quark states 



considered here are mass eigenstates, i.e. Cabbibo-Kobayashi-Maskawa mixing |12 
is suppressed for the sake of a compact notation. To simplify notation we suppress 
the obvious flavor indices V^^ — >■ in the following. 

A slight modification arises for the 7 and interference terms. 



a 



dE'd^l EM (Q2 + (M2o)2)g2 



m: 



')(!) 



j^py-M.7,^°)(2) 



int. 7,^0) 



/r fint. 7,Z0){l)u/(mt.7,^°)(l) , r (int. 7,Z0)(2) TT/(int. 7,Z0)(2)\ 

^ant.7.^")(i) = ^(fc^i (4^.7)(o))^ |A:V)(A:V|jf '^")(0)|A:a) , 

cr' 

^(int.7,z0)(2) ^ ^^^^1 (,(^.^")(o))^ |fcV)(A:V|jf -)(0)|A:a) , 

cr' 



r(mt.7,Z0) ^ r(int.7,Z0){l) , r (int. 7,Z0)(2) 



(9) 



Using Lorentz covariance and time-reversal invariance the nucleon scattering tensor 
W^f '"^ (g, P, S) can be expressed in terms of 14 structure functions (In our sign 



convention all e structures are multiplied by —i 



^q^q^\Fr\x. Q') + ip'q" + P^q^ ^^r\x, Q') 



V 
V 



—le 



PaS/3 {Pr), , T?{Pr) / 



V 



2v 



+ + ^1 ^ar\x,Q^)+p^p^^ar'^ 



q 



.Sq (Pr). ^2 



Wq''^af^\x,Q') + {p^q'^+fq^) ^,at\x,Q') 
+ [S^pT + S^r) ^€''\x. Q') + (S^P" + S^P") ^h^t\x, Q')- (10) 

Here u = p ■ q and x = /{2u). Furthermore the vectors p^i = Pfj, + qfit^/Q"^, and 
= + qfj, S ■ q/Q'^, which vanish if contracted with q^, are introduced. Note that 
terms proportional to q^ or g*^ give contributions of order mi/E when contracted 
with the lepton tensor L^^. The decomposition yields seven structure functions for 
parity-conserving processes and additional seven for parity-non-conserving processes. 
All 14 are given in table |. Six structure functions (-F4, -F5, (73, CI4, 05, ^2) are related 
to non-conserved currents and are non-zero only from non-vanishing quark masses. 

In the framework of OPE the moments of structure functions can be related to 
nucleon matrix elements of local operators with well defined twist. To this end the 
nucleon scattering tensor VF^f is related to the virtual forward Compton scattering 
amplitude T^f ''^ 

W^^^^ = ^(T^^^\qo + ^e)-T^^^\qo-^e)) , (11) 



Sq 



A-ni 



TL^'\qo± 



{2Ti)H\px-p-q) 
Pxo -Po- (go ± ^e) 



T.ipS\{jrf mX){X\jl^^\0)\pS) , (12) 

X 



iTii J UJ' 



n+2 



(13) 



with uj = 1/ X = 2z//Q^. In the following we deal with crossing-even and crossing-odd 
amplitudes 

T;Z\-q)=Til\q) , 

Tfr\-<1) = -T?r\<D ' (14) 

where T^^^"^) = T^Jf ") ± Tjyf^\ Note that the formal equation (|l|) is only valid for 
certain moments of the structure functions because its derivation requires for closing 
the integration contour according to fig. 2 for the n-th moment 

Tl£''\^) = -{-T Tir/\-u^) . (15) 

This condition if fulfilled in the case of the crossing-even amplitudes such as Tj^J , 
T^f°^ and T^^^'^) for the odd moments (even n) of the structure functions F2, F4, a2, 
04, &2, 91, 92 and F3 and for the even moments of Fi, F5, ai, 05, and 93. For the 
crossing-odd amplitude T^^""^^ the situation is just inverted. Furthermore performing 
the integration in equation ( |T3|) is possible only if T^i^^u) vanishes fast enough for 



u! —>■ 00, i.e, when the integration over the circle in infinity in figure 2 does not con- 
tribute. This is true for n > 1 but not necessarily for n = 0. (In practice one expects 
nonsinglet sum rules to converge better since there is no pomeron contribution and 
nonsinglet structure functions are less singular for small x.) 

To analyze the Compton forward scattering amplitude 

T^^'\q,p, s)=ij d'i^^Hps\T ((jr^)nojrHo)) \ps) m 

we insert the nucleon currents Eq.(R]) which we now write in the general form 



J/.//'(0 = g/(Or,//'g/'(0 • (17) 

Note that T^//' = T^jif. Restricting ourselves for the moment to the handbag 
diagram (cf. figure 3) we write 

+qf,iO)T,ffnSfiO,0^,ff,qjiO}\pS) • (18) 

Here Sf{^,0) is the fermionic propagator in the external field approximation for a 
quark with fiavor /. In the Schwinger formalism for calculations in external fields as 
proposed by Shuryak and Vainshtein 0, (p!8D can be written as 

9/(0) 



T,.uff' = -{pS\ |g/(o)r^//' ^ _^^^_^^ r.//' 



+g/-(0)r.// y _ _^ T,ff,qf,iO)^ \pS) , (19) 
with = id^ + gA^. Such an expression is very conveniently expanded in powers of 

i/Q 

1 1 °° / \\^ 



This is the starting point of our OPE analysis. Inserting the expansion in ([T9|) it 



immediately becomes clear that for / = /' all odd terms in j are symmetric under 
exchange of /x and u and thus contribute to the unpolarized structure functions, 
while all even terms in j are antisymmetric contributing to the polarized structure 
functions. 

Finally, to define twist the different Lorentz structures are classified according to 
their spin. Although this is standard in principle some ambiguities enter for more 
than two Lorentz indices. To avoid confusion we state explicitly the decomposition 
we use in appendix A. 
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III. PHOTON AND EXCHANGE 



In this section we perform the decomposition of the virtual forward Compton scat- 
tering amphtude for exchange. The corresponding expression for 7 exchange can 
simply be obtained by setting Vf ^ Qf and Aj ^ 0. For parity-even quantities 7 
and Z° exchange always interfere. The interference term is given by 

7-(int7,z") _ f7-(7+^o)A _ (r(^)\ _ (r(^'')\ (O-l) 

Af-.Af Af^O Af^Af 

We will now proceed as in and |15] except for using a different classification of 
matrix elements. The virtual forward Compton scattering amplitude for 7 and 
exchange reads 

+ imiAVf + An.U — ^ — -ri,[Vf + An,]qf{o]\pS) • (22) 



Expanding in — mf)/Q according to ( pOD we get for the zeroth order term 



T^r'^ = -E(P^I ^/(0)7.[^/ + ^/75];|74V/ + A/75]g/(0) 

0. order / 

-QfiOhAVf + Aj^,]^^,[Vf + An,]qf{0)^\pS) 
= -2ie,^J^ua%'\vf + A)) 

-2ie,,xJ^u:[2af\2VfAf) , (23) 

where the identity 

was used and the following matrix elements were defined: 

2a%'^S'' = {pS\qf{Q)^"^,qfmvS) , 

24°'V = (p5|g/0)7'^g/(0)b5) . (25) 
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The classification of matrix elements is to be understood as follows 



a 



im,t) . 
/5- • 



(26) 



m = number of — operators. 
t — twist. 
/ = quark flavor. 
5 = indicates 75. 

— = indicates odd chirality, left out 
in case of even chirality. 

The spin of the operator is s = 3 + m — t. The chirality of an operator is said to be 
even if 



(27) 



holds. The other case, with the positive sign, indicates chiral odd operators. The 
contribution of first order is symmetric in /i and u. It reads: 



1. order 



f ^ g g 

+?/(0)7.[^/ + An^]iir - mf)^j,[Vf + Af^,]qf{0)^\pS) . (28) 



q^ 

We decompose the above expression into operators of well defined spin, which are 
parametrized in terms of reduced matrix elements. We write the result in a way 
matching the decomposition structure of the nucleon scattering tensor Wfj,,^. In this 
way we get 

M2 



7^(7,^°) 

1X1/ 



= -Qnu H — 

1. order \ q 



1 1 



, / , Sq 



Sq 

2VfAf + q^q,-^ 



M . 



2VfAf 
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Sq 



+ ^iS^Pu + SuPf,) 
,(0,2) 



2VfA 



. 1,2 , . rUf 0,2 



2VfA 



2VfA 



(29) 



For the matrix element af^_ we employed the equation of motion. The matrix 
elements are again classified according to Eq.(p6D. The exact form of the operators 
can be found in appendix ^. 

The contribution of second order in 1/Q is antisymmetric in fi and u. It reads 



= -E(P^I^/(0)7M[^/ + ^/75];|(f -rn/);|(f -m;)-|74V/ + ^/75]g/(0) 

2. order ^ q q q 

-M^hAVf + ^/75]^(f - m;)^(f - mf)^^^,[Vj + A^75]g/(0)b^) 

(30) 



and can be decomposed into 



2. order 



-- E -8^e^.A.gag/3g'(p5|g/0)P-PVg/(0)75b5)(^/ + 4) 
+VgA,^^^^,<x^^^(p5|^^(0)p-ppyg^(0)|p5)(lf + 4) 
+2g2^^g^e^,A.e%5(p5|g/(0)7Vi^'g/(0)b5)(K; + 4) 
+Vmyg„(p5|g/0)P"[7;.7. - lul,]qfmpS)A} 
-Ate,,,^qVm}{pS\qf{0)ri,qfmpS)iV} + A}) 
-8ze^,Mpq\pS\qf{0)P"pfrqfmpS)2AfVf 



+4qVe,uX.e\^^{pS\q^{0)P^P^^^^,\pS)2AfVf 



(31) 



Here we use the following relations that result from the equation of motion ^qf 

Vfifqf 

m;(p5|g/0)7^g/(0)bS) = (p5|g/0)P^gX0)b5) , 
(p^|g/0)75g/(0)|p5)=0 , 
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{pS\qf{Q)P^^,qfiO)\pS) = , 
(p5|g/0)(P^7.7A + + ^A7^7.)75g/(0)b5) = . 



(32) 



The totally antisymmetric spin-0 part has to be treated separately. It can be trans- 
formed in the following way (Note, that these equations are only valid when sand- 
wiched between {pS\qj{0) and qf{0)\pS).): 

e''^''^P^Pp^, = i{P^-m))r^, , (33) 
e"^^'<P^Pp^,^^ = i{P^-m))Y , (34) 

e""^^7a7/3P7 = 2^^/7''75 , (35) 
e""^^7a7/3i^775 = . (36) 

With the definition of the dual gluonic field strength tensor igG^r = "^f ^o-tq/jG"^ = 
\^uTai3[Pai Pp] and the equation of motion the operators P^7'^ and P^7'^75 can be 
transformed into a mass component and a gluonic component. In terms of matrix 
elements this transformation reads: 



V5 



(37) 
(38) 



where the reduced matrix elements are listed in appendix ^ Here also a second 
relation is derived: 



(2,3) _ ~(2,3) l^f (1,2) 



V5 



V5 



(39) 



Collecting all terms we get for the second order of our expansion in {f — mf)/Q: 



q^S'' 



-It 



2. order 



3 -'f' 
8 1 
'3Q 



^ 3 (2,2) , 4 3^(2,3) "ilVl- / (2,2) 

3 



4M2 / 
V5 +q7^^l 



(0,2) 2 
2^«/5 



9Q^ 



. - , I /l~(2,3) , A~(2,A) 
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4 m. 



'f, ,„(0,2) 



{Vf-A)) 



q p 



qXpCr 
2^ 



/ ' q r)2 V J J 



9 Q 



9Q' 

(0,2) 2 



~(2,4) 



2VfA 



(40) 



IV. VF=t-EXCHANGE 



In the case of W exchange we use the current 



(41) 



where the eUipsis denotes currents involving heavier quarks. For simphcity we will 
discuss only the m — c? current. For other flavors the relations look just the same. 
The equation analogous to (^) is then given by 

T^r^ = ^ I d'^e''^HpS\T(u{Ol,iV + A^,)d{Omiu{V + Aj,)u{0)) \pS) 



-{pS\^u{0)^^{V + A^,] 
+d(0)7.(V^ + A75)- 



1 



1 



^^iV + A^,)diO)\\pS) 



(42) 



where the expansion in If —mf becomes more complicated since two different flavors 
are involved. In each order symmetric and antisymmetric terms appear. For the 
zeroth order term of the expansion we have 



0. order 



-{pS\{u{Oh,{V + A75)^7.(^ + ^75)^(0) 
+miu{V + A^,)^^,{V + A^,)dm\pS) 
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-9iiu + 



Lua. 



(0,2)1 



2^(0,2) 



-Loa, 



(0,2) 
m5 



2u 



u) a. 



V 



u 



+ 



qp^qv\ Sq (0,2) , 1 



V 



■If] 



2v 



[2ai°'^)] I + (m ^ d, ^ z/, g ^ -g) 



(43) 



The fact that in first order of the expansion we get contributions to current non- 
conserved structure functions, namely ^4,^5, 62, and is due to an incomplete 
definition of T^^. Indeed, in zeroth order no quark masses appear and conse- 
quently, both the vector and the axial vector current should be conserved, requiring 
q^Tiiv\'i. order = 0. The tcrms proportional to F4, F5, 62, and 04 are therefore unphys- 
ical and must be subtracted. In this sense our result (|43|) has to be rewritten in a 
form where only current-conserved contributions occur: 

q^S"" (0,2) 



rpiW-) 



0. order 



-9i,v + 



qt^qv 



q 



'^^fiuXa ' 



-Lua 



u5 



{V + A' 



'9^cu 



2u 



-te^uxa^uj[2a'^'^'>] ^2VA+{u^d,fi^u,q 



-q) 



(44) 



This phenomenon arises only in the first order of expansion. Formally it comes 
from unphysical contact terms (seagull terms) generated by derivatives acting on the 
theta function in the time ordered product [0. In this way we obtain, expanding 
the propagator to first order 



m(vy-) 

^1/ 



= {pS\{mi,iy + ^75)^(f - m,)^7.(V^ + ^75)w(0) 

1. order q q 



+rf(0)7,(V + A75)^(f - ma)\^,{V + A^,)d{mpS) 
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(0,3) 



1 



1 om 



3'"w^(0,2) 
u5- 



--tf — — a 
2 M 



u (0,2) 1 



M 



u5- 



2VA 



O^. (1,2) _ (0,2). 

2 V^us ^ "us- / 



-^6,...^[2a;^a^^)]2yA 



M 



(45) 



The expression for the second order term 



2. order 



+d(0)7,(y + A75)^(f -m„) 



,7M(y + >l75)^^(0)}b5> 



is again decomposed into components antisymmetric and symmetric in /i and u. The 
latter one reads: 



rp{W-) 



2. order 
sym. in (xi/ 



3^(2,2) 8 M\ r (2,2) ^-(2,3) ~(2,4)i 



a; a, 
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+ 



4muM, ,„(l-2) 
3 Q2 ■ 



(0,2) 



m^, ,3^(1>2) 
u5- 



M 



uj a. 



1 m„, ,4^(1,2) 



2 M 



a; a. 



«5- 



2 Q2 ["^w 



21, ,2^(0,2)1 



m J a; a 



u5 



2, ,2^(2.2) 



^6j2o(2,3) , 2HhLa,2„(1.2) 



+ 2i;('^mP- + '^-P/^) 



4, ,2„(2,2) , 8, ,2~(2,3) , 8 m^M^ (1,2) 



u5 



+ 



(0,2)' 



3 



(1,2) 4 1 r 2 _ Q™21_IU, 
'■u5- ^ 3Q^[''''u •^'''■dJ"«5 



+ + [-2^(^2^(1,2) _ |n^«(i,2) _|^[^2 _ Sm^aS')] } 

+ {u d, /J, u, q —q) 



(46) 



The corresponding antisymmetric part differs from the case only by the mass 
terms and, of course, by a factor 1/2. 

q^S" rl . (9. 9.^ 2 



rp{W-) 



2. order I jy 

asymm. in ni^ V. 



3 (2,2) , ^ 3~(2,3) 



,r«^2'^) ^ 4;i(2,3) 4~(2,4)n 1 1 r 2 , q^2l^(0,2) 



q^p" 



^ 3r (2,2) ~(2,3)i 



■^^[26m^ + 18m^]ai°'2)u; 



2y>i 



+{u <-> d, // <-> i/, g <-> — g) . (47) 

We have discussed matrix elements up to twist 4. At this twist level the cat ear 
diagram (see figure 4) may contribute. It describes non-perturbative interactions 
between the hit quark and the spectator quarks in the final state. However, it 
turns out (see appendix B) that the matrix elements of the cat ear diagram yield 
contributions of structure functions that are one order in higher such that we 
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do not have to include them in the order we are working. Also, due to symmetry 
properties cat ear diagrams can only contribute to Fi, F2, ai, and 62, as discussed in 
appendix B. 

V. SUM RULES 

Having analyzed the virtual forward Compton scattering amplitude completely up 
to order [(^ —mf)/Q]'^, we have determined the lowest moments of the 14 structure 
functions. The results are compiled in table II and the resulting sum rules are listed 
in table III on the leading twist level. The mass corrections to these sum rules can 
be read off from table II so we do not repeat them. Table IV gives some relations 
between structure functions, valid in the limit — > 00, i.e. neglecting higher-twist 
corrections. In general, here we do not consider possible anomalous contributions 
arising from the divergence of pseudovector currents. 

Let us add some comments to these tables: 

1.) The index / runs over all quark and antiquark flavors, i.e. f — u,u,d,d, . . . 
Due to the opposite helicity of quarks and antiquarks we have 

Vf^Vj ; Af^-Aj . (48) 

Furthermore, we use in the case of neutrino scattering for the weak coupling constants 
the abbreviation 

and for brevity, in the case of - exchange, we introduce the following flavor 
combinations: S (singlet), V (valence) and furthermore the combinations and 
AV. 
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as ■= Vwiia-u + a-u) + {cid + a^) + (further generations)] , 
av '■= Vw[{(iu — ciu) + {cLd — CLj) + (further generations)] , 
OAS := ilw[{.<iu + O'u) — {cid + 0-^) + (further generations)] , 
o-Av '■= Vw[{(iu — ciu) — {cLd — 0-^) + (further generations)] . (50) 

The relations between these flavor combinations and the symmetries of the ampli- 
tudes are presented in table [V]. 

To simplify the mass terms we define 

{nis - m's) as = [(m„ - m^) a„ + {niu - m^) 

+ (^d ~ "^u) ad + {rrij — m-u) a-^ + (further generations)] , (51) 

and analogous expressions for V, A S, and A V. 

The constant C appearing in table III for the polarized Bjorken sum rule has the 
value 

C= ^ +Tr-^[4sin^;^-l] + |^^[2sine^-l] . (52) 
3 3 sm 20w 3 sm 29w 

^ ' ^ ' 

"f—Z^^ interference 

Isospin T3 and baryon number B are given by 

E {af'^-af'^)- E (ar-af))=2T3 , (53) 

f=u,c,t f=d,s,b 

E (ar'-«r)+ Y. («r-4°'')=3S , (54) 

f=U,C,t f=d,S,b 

Finally we note that we do not include here the radiative corrections and chose 
for simplicity the normalization point to be /i^ = Q^. 



2.) The twist-2 matrix element of the tensor operator {—ija^i, determines the 
transversity distribution hi |T7] and does not appear at leading twist in deep in- 
elastic lepton nucleon scattering. As a chiral-odd operator it is related to spin-flip 



18 



processes which are always suppressed by the quark mass. As may be expected be- 
forehand this operator determines the polarized structure function which is related 
to a non-conserved current and vanishes identically for massless quarks. 



3.) The simplest scalar twist-3 operator is the scalar operator 



{pS\qf{0)qfmpS) = 2afl'^M 



(55) 



which for the flavor combination ^{mu + mii){uu + dd) is known as the cr-term, giving 
rise to a sum rule for F4 which, however is expected to be divergent due to a Regge 
analysis Since the scalar operator is chiral odd it will appear in deep inelastic 
scattering only if quark masses are not ignored. 
4.) The twist-3 quark-gluon-quark correlator 



{pS\qf{0)- 



c 



qf{Q)\pS) — traces 



2a 



(2,3) 
/5 



-{2paPf3S„ - PpPaSa - PaPaSp) - traCCS 



(56) 



(often referred to as S'^^ ) contributes to the first moment of gi and the third moment 
of g2. This operator contributes also higher- twist terms to a\ and 02- 



5.) The twist-4 quark-gluon-quark correlators 

{pS\qf{Q)gG"^^pqfmvS) = 2a%'^M^S'^ . (57) 

contributes to gi and ai and this matrix element is commonly referred to as /*^^^. 
The spin independent twist-4 operator defined as 

(p5|g/0)(7G"^^75g/(0)b^) = 2af^)MV , (58) 

gives a higher-twist corrections to the Gross-Llewellyn Smith sum rule. (The matrix 
element of the analogous twist-3 operator vanishes due to its symmetry properties, 
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i.e. af^'^ = 0.) 



6.) It is possible to relate the higher-twist contributions for the second moment 
of g2 to the structure function and to obtain an interesting relation, which reads 



xdxg 



{u-u) 



2 Jo 



{u-u) 



(59) 



/o - 2 JO 

7.) We cannot confirm the relation a2(x) = 2xai{x) for —>■ oo, which was 
claimed in but rather reproduce the relation found by Ravishankar [|l8l and 
earlier by Dicus [|T9[] : 



2xai(x) = a2(.x) + + b2(x) 



(60) 



8.) For the polarized parity-non-conserving structure functions we reproduce the 
relation derived by Wray [pO] by means of the light-cone quark algebra: 



dx 



-2%^ + 0{1/Q') 
9v 



(61) 



While the current non-conserved structure functions, i.e. 04, 05, and 62, will hardly 
be measurable in practical experiments, the current conserved structure functions 
ai, 02 and 61 should be experimentally accessible 0. They give interesting additional 
information on the structure of the nucleon. 



9.) We would like to stress that the operators connected with 0^5^'' and 



,(2,4) 
^5 



measure important properties of the nucleon, namely contributions of the collective 
color electromagnetic field to the spin and to the momentum of the nucleon. In the 
rest system = (0, S), = (M, 0) the spin-dependent operators can be written as 
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{pS\l-B\pS) = , 

{pS\ [-Ba£ + (ja X Ea)] \pS) = 2M'af^^S 
{pS\ \2Bafa + il X E,)] \pS) = m^a%'^S 



(62) 
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where = —gq'j^taq denotes the quark current, B'^ the color magnetic field, and 
the color electric field, af^'^^ and were calculated in relativistic quark mod- 



els such as the MIT bag model as well as in the framework of QCD sum rules 
2^,|23|pl| making it possible to estimate the strength of the color electric and mag- 



netic fields within the nucleon separately. 

Analogously, for the axial vector current j^g = —gqj^j5taq we get 



{pS\l, ■ B\pS) = - 
{pS\ + (^5 X Ea)] \pS) = 

{pS\ [25,J°5 + (^5 X Ea)] \pS) = 



(63) 



The two operators jas x Ea and Baj^ determine a vector direction. As in the nucleon 
rest frame the only distinct direction is given by the spin, which is an axial vector, the 
matrix elements of ja5 x Ea and -BaJas niust vanish. We would like to point out, that 
the only nonvanishing matrix element in this context, namely CLf''^\ determines the 
higher-twist corrections to the Gross-Llewellyn Smith sum rule and was calculated 
in the QCD sum rule approach in [EH . 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have presented a complete analysis for the first moments of the 
structure functions of the nucleon scattering tensor for 7, Z^, and exchange. All 
our calculations include the effects of finite quark masses. 
Besides reproducing well-known results we found the relation 

62(3^) = 2xa5{x) = 4x^a4(x) . (64) 

The extension of our work to the third order in ^ — n^'f will require a consistent 
treatment of the cat ear diagrams. 
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APPENDIX A: DEFINITION OF THE MATRIX ELEMENTS 



For the sake of brevity we introduce the notation of the totally symmetric tensor 



I ^cti 



I.e. "Sqj^c 



vanishes for the same index combinations as 



and whenever there is an index combination such that e^jc 



±1, then 



1. All other conventions are taken from \14 



i) g'^'^ has spin 0. 

ii) Each Lorentz vector and 7^^ has spin 1. 

iii) The operator can be decomposed into a symmetric and traceless spin-2 
part, an antisymmetric spin-l part, and a spin-0 part according to 



1 



J Spin- 



+ 











+ 





• (Al) 



iv) The operator P^P^^'*' is decomposed into spin-3, spin-2, spin-l and spin-0 
parts. There is one ambiguity in the definition of the spin-0 part: Sandwiched 
between nucleon states, it can be transformed into a spin-l operator via the 
equation of motion (see Eqns. ( p3H35|) ). 



papP^ — ["p"p/3^7 



7' 



Spin— 3 



+ 



pO'pP^ 



Spin-2 



+ 



Spin-l 



+ 



papP^' 



Spin-0 
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6 



18 ^ 



- ^2P" P^7''' — P'^ P"7^ — P^ P'^j°''j 

-- (5"'3[2P27T -fP^- P^f] + 5"^[2fp/^ - P'^f - P^j^ 
9 ^ 

+^7/3 [2P"^ - f P" - P^7"]) 



+ 



+ 



— (^"^[5P^7^ - f P^ - pTf ] + t,"T[5f P^ - P^f - P27^] 
18 ^ 

+^7/3 [5p«^ - f P" - P27°^]) 



^e"'^^^e5,,,P<^P^7' 
6 



(A2) 



v) For the operator j°'j^P'^ the situation is essentially the same as in iv), however, 
this operator has no spin-3 part due to the fact that different gamma matrices 
anticommute. 



a^P p7 ^ [7^7'^ 



7 7' 



Spin-2 



+ 



Spin-1 



+ 



"^/3p7 



7 7 



Spin-0 



- (27"7^P^ - 7'^7TP" - 7T7"P'^ - 27^7"P^ + ^^^f^P"" + 7"7Tp^) 
7^" - 7''f ] - 7" - 7°f ]) 



+ 
+ 



(A3) 



According to this decomposition we define the various matrix elements in the 
following manner, ordered according to the number of vector indices. 
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Scalar and pseudoscalar operators 



{pS\qfiO)qfmpS)^2afl'^M , 
{pS\qf{Q)^,qf{Q)\pS) = . 



(A4) 
(A5) 



Vector and pseudovector operators 

{pS\qf{0)P,qfmpS) = mf{pS\qf{0)^,qfmpS) , 
{pS\qf{0)P,^,qf{0)\pS)^0 . 

The last equation was simplified using the equation of motion. 

{pS\qf{0)^,qfmpS)^2af^p, , 
{pS\qf{Oh,^,qf{0)\pS) = 2af^^S, . 



(A6) 
(A7) 



(A8) 
(A9) 



Rank-two tensor operators 

In the spin-2 case the parametrization of the two possible operators is straightforward: 



,(1,2) 



{pS\qf{0){P,^, + Pul,)l,qfmpS) = 2^^ ) \p,S, + P.S, 



, (AlO) 
(All) 



For spin-1 only the pseudo-tensor operator yields a non-zero matrix element: 

(p5|g^(0)(P^7. - P.1MQ)\PS) = , (A12) 
{pS\qf{0){Pi,7u - Pvlix)lbqf{^)\pS) = mf{pS\qf{0){-ia^^)-f5qf{0)\pS) 



(A13) 



af^f! is connected to the first moment of the parton distribution hi{x) which is 



proportional to the structure function Gs{x). 
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Rank-three tensor operators 

The operators PaP/sJa and 'jalpPa have to be considered. 
lalpPa- Spin-2 

Again, only the operator with 75 has a nonvanishing matrix element, i.e. 



{pS\ qf{0he.7f}Pa ?/(0) \pS) = , 



(A14) 



{pS\qf{0)jalpP, . _ 75g/(0)|p5) 



{pS\qf{0) 



_ ^ (1,2) 



spin— 2 



(A15) 



aW is connected to hi{x) as indicated by {—i)(Tai3- 



lalpPa-- Spin-0 

By means of the equation of motion we get: 



3 

2 



^(0,2) 



-{pS\qf{0)e^p,se'rpxYYP'lBqfmpS) = . 



(A16) 
(A17) 



PaPp^a-- Spin-3 

The parametrization of the spin-3 part is straightforward: 



{pS\ qf(0)P^Pp^, 



2a 



(2,2) 



PaPpPa ^[gapPa + QpaPa + QaaPfi] 



(A18) 



{pS\ qf{0)PaPpJ, spin-3 ^5^/^°) l^-^^ 

25 



2a 



(2,2) 
/5 



(A19) 



PaPp^a- Spm-2 

In the spin-2 contribution the gluonic field strength tensor comes into play due to 
the following two identities: 

^(p/3pa_pap/3^)|^^ , (A20) 



and 



I ^2P° P^'y'^ + 2P^ P'^7'^ — P^ P^7" — P'^ p^7" — p'^p°'^l^ — p^'p" 



.1 p2 

3 



— m 



- (27^7^^ P'' — 7'^7° P^ — 7°^7^ P'^ — 27'^7^ P"' + 7^7'^ P^ + 7^7°^ P'^^ 



^27^7"^ P^ — 7"^ 7^ P° — 7^7"^ P'^ — 27'^7'^ P^ + 7^7"^ P" + 7^7^ P'^^ 
-^(/"[f7^-7'^f]-/'^[f7"-7"f])l|75 . (A21) 



Again, one of the two possible matrix elements is zero 



(^S\qf{Q)P^Pp^, ^^^_^qfmpS) 
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{pS\qf{% 
1 



+ - 







'^gapgGaSl^ - gpagGaSl^ - gaagGpSl^ 



75 



(A22) 



This implies that the gluonic matrix element is zero as well. In the other case we 
have for the gluonic representation and for the representation in terms of covariant 
derivatives: 



1 



6 



1 

+3 



2a 



{lagGpa + iugG 
'^gangGasi^ - gpagGasi^ - gaagGfssi 

(2,3) 



/5 



1 

-{2paPf3S^ - PpPaSa-PaPaSp) —[2gaf}S^ - g^aSa " gaaSfs] 



, (A23) 



and 



{pS\qfiO)P^P^'y„ . 759/(0) 



2^(2,3) 



spin— 2 
1 

■^{'^PocPpSa - PpPcSa - PaPaSf)) —[2g^pS^ - QfjaSa " goc<jSp] 



(A24) 



These two representations are equivalent except for mass terms. Explicitly one has 
the identity 



(2,3)_~(2,3) rUf a^2) 
«/5 - «/5 - 'J[^f5- 



(A25) 



Spin-0 

For the spin-0 contribution there exists also a transformation between a gluonic rep- 
resentation and a representation in form of covariant derivatives only: 

27 



2ie 



afiaS 



V5 



2„(0,2) 



, (A26) 



and analogously 

= ieapas{pS\qfmP' - m})j'qf{0)\pS) 



= 2ie 



Obviously 



aj3a5 



,^2 (2,4) ,^2~(2,4) 2 (0,2) 

,^2 (2,4) ,,2~(2,4) 2 (0,2) 

M ay ' — M ay ' — irifaf ' 



2 (0,2) 



(A27) 

(A28) 
(A29) 



APPENDIX B: OPERATORS GENERATED BY THE CAT EAR 

DIAGRAM 

In the required approximation the cat ear diagram (fig. 4) is not connected with 
the external gluonic field, so it is sufficient to analyze ordinary Feynman diagrams 
in zeroth order in — mf)/Q. In the case of 7 or exchange the corresponding 
virtual forward Compton scattering amplitude then is given by 

= T.t/ I / dz,dz2(j)S\{qyi)^^{Vf + Af^,)isyi,z,)Yt%{z^) 

rder ff ^- J J ^ 

+qyz,)YtSSyz,, 07m(^/ + ^/75)g/(0}^^S(^i, ^2) 
X |g^,(^2)7"^'^<5/'('^2, 0)7.(T^/' + ^/'75)g/'(0) 
+g/'(0)7.(^/' + ^/'75)^^/'(0,^2)7'^^V(^2)}|p^) , (Bl) 
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rpicaX ear 7,2") 



where we have to insert the free massless fermion-propagator. D^{zi,Z2) is the 



gluonic propagator connecting the "ears" of the cat ear diagram, ta, a 
denote the color matrices (tr(t'^t'^) = 5^^/2). In that way we obtain 

= -yE(P'5|k(0)7M(V} + ^/75)j7pt'^g/(0) 

0. order ff ^ ^ 



T-i(cat ear 



1 



(0) 



~9 



crp 



X 



?/'(0)7.(^/' + ^/'75)^7.^V(0)}|P^) 



,8 



(B2) 



The corresponding expression for the W exchange is 



^^(cat ear W ) 



0. order 



-u{0)jpt''-rj^{V + Aj,)d{0) 



X rf(0)7,t"-7,(l^ + A75)«(0) 



-c^(0)7.(^ + ^75)^7.*"«(0) \\pS) 
+ further flavor combinations 



(B3) 



(The case of exchange is obtained by exchanging the quark flavors.) Introducing 
reduced matrix elements of the occurring four-quark operators we finally get in zeroth 
order 



T-i(cat oar 



^9^1 \ -9nu + '^^ 

0. order I V q , 



+ 



Sq 



-9ny 



q,,qv 



,2 



UJ^ ^(2) 



V 



2— 



2—C 



(2) 
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a-T(cat ear W 



0. order 



Q2 «rf=^ Q4 



^(0) , ^(2) 



1/ 



9 ^(2) 



9 ^ r^(2) 



(B4) 



where only the expression for the u d flavor combination is given. Here the flavor 
independent matrix elements C^j^ are defined by 



Mi) \ " 



=/x •- Zl^^//=/x' 

/ 



J = 0,l,2,... . 



(B5) 



Note, that the cat ear diagram (|B4| ) is symmetric in fi and u. Therefore in zeroth 
order oi If —irif it gives only twist-4 and twist-6 contributions to structure functions 
associated with symmetric Lorentz structures, such as Fi and F2. Since in the re- 
quired approximation quark masses do not occur, the cat ear does not contribute to 
structure functions given by non-conserved currents. Therefore we find that cat-ear 
diagrams contribute only to the structure functions Fi,F2,ai, and 62- 

For the definition of the reduced matrix elements of the cat ear diagram we need 
the following definition of four-quark operators: 



^xxjf •- 



qf{0)ri,t%iO)qj,{0)r'l5t%{0)AfAf, , 
qf{0)ri,t%iO)qf,iOh'''t%{0)AfVf, , 
qfiO)rt^qfiO)qf,{0)r'l5t%{0)VfAf, , 
qjiO)Yt%iO)qj,iO)r't''qfmfVf' • 



(B6) 
(B7) 
(B8) 
(B9) 



Note that for 7 and exchange / = /', and in the case of exchange we always 
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have Vf — Vf — V and Af — Af — A. We then define the reduced matrix elements 
to be 

-J^^^'w(055j/' + ^xl//')|P^) , (BIO) 



= (:pS\\{6gjf, + 6^;,^^, - - <^xxjf') \pS) , (Bii) 



2M^cgL = (p^l w(055j/' + ^xl//') \PS) , (B12) 



-]g''^'9aa'id^^,fr + dTifr)\pS) , (bi3) 



= {pSlkdrxjr + OTifr - O^^hr - ^xhr) , (bi4) 



0^{pS\g^AO^^',ff'+OTlff')\pS) . (B15) 

Note, that the operators Oxxjf',05Xjf',Ox5jf' can be transformed by means of 
the equation of motion 

D'^Gl^^-Y.gqf^^tyf, (B16) 
/ 

where Da is the covariant QCD derivative and the gluonic field tensor. However, 
as this operation cannot be applied to C*55,//', we did not perform it to keep a uniform 
notation. 
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Figure Captions 



Figure 1: Kinematic variables in deep inelastic scattering: An electron with four- 
momentum k and spin a is scattered by a nucleon with four-momentum p and polar- 
ization A through the exchange of a gauge boson with four-momentum q. The nucleon 
fragments into a variety of hadrons X, which together carry the four-momentum px- 

Figure 2: The contour of integration of the virtual forward Compton scattering 
amplitude T^jy. The regions where the function T^i,{cu) has singularities are shad- 
owed. The four big arrows show the integration path. It can be completed by the 
outer circle in the infinity thus being equivalent to the inner circle. 

Figure 3: The hand bag diagram: It represents the virtual Compton forward scat- 
tering amplitude, for quark (solid line) gauge boson (sinuous line) scattering. 

Figure 4: The cat ear diagram: The lines represent: Gauge boson (sinuous fine), 
gluon (curly line), quark (straight line). 
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TABLES 





unpolarized structure functions 


polarized structure functions 


parity-conserved: 






current conserved 


-^15 F2 


91, 92 


current non- conserved 




93 


parity- violating: 






current conserved 


F3 


Ol, 02, h 


current non-conserved 




04, 05, b2 



TABLE I. The 14 nucleon structure functions. 
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Moments of structure functions 



(7,2") 



(7,2°) 



J^xdxF^ 



(7,2°) 



(0,2) 2 M2 
9 Q5- 



Y.faf^hVfAf + 0{l/Q') 



Ef^aff^Aj + Oil/Q') 



+ 0(1/Q2 



+ 0(1/Q2) 



J.'dxFf'^ 

xdxFi""^"^ 
lo'x'dxFt''^ 



xdxF^ 



lo'dxFt"^ 



J.'xdxFf^ 



x,Q^) = 



ag^)+0(l/Q2) 



= 2ag^^+0(l/Q2) 



,(0,2) 4M2ro„(2,2) 



(2,4) 1 



-2a(,^|)+0(l/Q2) 
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Jx'^dxFf^''\x,Q^) = + C>(l/g2) 

/x3dxFf-")(x,Q2) = + 0(1/Q2) 

/o^ dxFf -'^^ (x, Q2) = +2^ - m'ly]a'-l^^ + 0{1/Q^) 

J^xdxFt^''\x,Q^) = + O(l/Q2) 

J^x^dxFt''\x,Q') = + 0(1/Q2) 



J,\^dxg^{''^'\x,Q') 



V + 1*^' r„(2.2) , . (2,3) . (2,4) .m, (1,2)\ 



E/Hf^(F/ + 4) + o(i/g^ 



+ 0(1/Q4) 



+ C'(l/Q^ 



= + C'(1/Q2) 



llxdxgi ''\x,Q'^) 
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Jq x'^dxg[''^''\x,Q 



ag'^+Oil/Q'] 



J^dxgt''\x,Q') 

j'^ xdxg^2 ''\x,Q'^) 
!^xHxgf^''\x,Q^) 



= + C'(1/Q4) 



,(1,2) 



,(0,2) 



2 '^AVS ~M^"Ay5- 



+ 0(1/Q2 



,(2,2) „(2,3) ms„(l,2) 



^55 



*55 



+ 0(1/Q2) 



j^dxgf-^'Xx.Q^) 
J^xdxgS'^''\x,Q^) 



+ 0(1/Q2) 



J^dx4''^°\x,Q^) = + 0(1/Q2) 
J,' x'dx4''''\x, Q2) = - k^affjAfVf + 0(l/g2) 

dxb^^''°^ {x, Q^) = -2 E/ ^4° - ^/^/ + ^( VQ^) 



dxa^^-'^Vx _ _a(0,2) _ 8 M2 r (2,2) _ 2 (2,3) (2,4), 4 mA^M (1,2) 

J^xdxaf^''\x,Q^) = -aS^f 
Jo^ x^dxaf-''\x, g2) = -ag|) + 0(1/ Q^) 



J,'dxar''\x,Q^) = + O(l/g2) 
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xdxa2 ^\x,Q 



f^x^dxaf+''\x,Q^) 
J^x'dxaf-''\x,Q^) 



J,Uxat"\x,Q^) 
xdxa'^^^''\x,Q'^) 

J,'x^dxat''\x,Q') 



+ 0(1/Q4) 



!Uxhf^^\x,Q^) 
Joxdxbf~''\x,Q^) 



4„(2.2) , 8_(2,3) 8mAS„(1.2) 



+ 0(1/Q2) 



Jq xdxb'^ ^\x,Q 



2^ae+0(l/Q2) 

«Akl + 0( VQ') 



TABLE IL Moments of structure functions with calculated corrections. Both the 
neutrino and the ^/Z^ case are given. 
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sum rule 


reference 


unpolarized 


dx{Fr - = IVW 


Bjorken (1967) 






Adler (1966) ^ 






Gross, Llewellyn-Smith (1969) § 


polarized 


Sldx{g{{x)-gUx)) = %'^ 


Bjorken (1966) @ 




fl ^^^P(^\ - 1 .9-4 9F/D-1 


Ellis, Jaffe (1974) H 






Wray (1972) Q 



TABLE III. Sum rules in the Bjorken limit (Q^ co). 
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relation 


reference 


unpolarized 


F2{x) = 2xFi{x) 


Callan, Gross (1969) [H 


polarized 


twist — 3 contributions 

+ 

and mass corrections 


Wandzura,Wilczek (1977) (H 




b2ix) = 2x05(3;) = 4x^04(2;) 






2xai(x) = 02 (x) + 61 (x) + 62 (x) 


Dicus (1972) |1| 



TABLE IV. Relations between structure functions in the Bjorken limit (Q^ 00). 





parity-conserved 


parity- violating 


amplitude v + 


(singlet) as 


(valence) ay 


amplitude V — v 


aAV 


OAS 



TABLE V. The four possible flavor combinations for neutrino scattering. 
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p, X 

Figure 1 




Figure 2 



Figure 




Figure 4 



